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Abstract 

We propose a formulation of the QCD partition function which 
leads to a series expansion for the quark determinant in any given 
baryonic sector. The r-th term gives the gauge-invariant contribution 
of the valence quarks plus r quark-antiquark pairs. This expansion can 
be used to investigate any baryonic sector, starting from the nucleon 
up to high baryonic densities. 
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1 Introduction 



A clearcut separation of the valence and sea quark contributions to the QCD 
partition function might help in the study of a number of issues, starting with 
the spin content of the nucleon |L] up to QCD at finite baryonic density 0. 
The identification of the different contributions can be achieved by means of 
an expansion of the partition function with respect to the number of quark- 
antiquark pairs, which is the goal of the present paper. But the practical 
usefulness of such an expansion depends on the number of terms necessary 
to get the desired accuracy 

In the case of the nucleon, for instance, only a few quarks of the sea are 
expected to contribute significantly. If the partition function is split in the 
corresponding terms only a few will be important. 

At first sight the case of QCD at nonzero baryon density appears quite 
different, because in principle one is interested in the termodynamic limit. 
To investigate this limit an infinite number of terms would obviously be nec- 
essary, making our expansion of little use. In actual numerical simulations, 
however, the physical volumes accessible are of the order of a few fm 3 . If we 
consider that normal nuclear density is about 3 quarks /6 fm 3 , with 6 valence 
quarks we are already at about twice the ordinary nuclear density, and it is 
reasonable to expect that the number of important quark-antiquark pairs be 
comparable to the number of valence quarks at least at moderate density and 
temperature. 

To our knowledge there exists no clearcut separation of the contribu- 
tions we are discussing. To tackle this problem we need, to start with, a 
formulation of the partition function in a given baryonic sector. One such 
formulation has been constructed by introducing an imaginary chemical po- 
tential in the grand canonical partition function and by projecting in a given 
baryonic sector by a Fourier transformation 0,0. But the result is not very 
transparent with respect to the distinction of the different contributions we 
are interested in. 

We consider here a different form of projection, which leads to a series 
expansion for the quark determinant. The terms of this expansion have the 
desired physical meaning: The r-th term gives the gauge-invariant contri- 
bution of the valence quarks plus r quark-antiquark pairs. In view of the 
above considerations, we hope that our expansion will have a reasonable 
convergence in present numerical simulations. 

The technique we employ is fairly general, requiring only the knowledge 
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of the transfer matrix in the Fock space of the quarks. But in this paper we 
restrict ourselves to the case of Wilson fermions. We hope to be able to treat 
the Kogut-Susskind fermions in a future work. 

We outline the derivation of the series expansion in Sec. 2 focussing on 
the strategy, leaving out many details of the calculations for Sec. 3. In Sec. 4 
we give our conclusions. 

2 The series expansion of the quark determi- 
nant in a given baryonic sector 

In this Section we outline the steps leading to a series expansion of the quark 
determinant in a given baryonic sector. As already said the presentation will 
be somewhat schematic. 

The quark determinant in the absence of any condition on the baryon 
number is 

detQ = f [dipip] exp Sp. (1) 

Sf is the quark action and Q the standard quark matrix which will be spelled 
out later. Our strategy is to write det Q as the trace of the transfer matrix 
acting in the quark Fock space, impose the restriction to a given baryonic 
sector, and then rewrite the trace as the determinant of a modified quark 
matrix. The round trip is done by mapping the Grassmann algebra generated 
by the quark fields into the Fock space following the construction of Luscher 
||. But while his paper is based on the mere existence of the map, for us it 
is essential a concrete realization by means of coherent states. 

In the transfer matrix formalism the (euclidean) time is treated differ- 
ently from the spatial coordinates. So we must also treat it differently, and 
we assume a convention of summation over intrinsic quantum numbers and 
spatial sites only. These have coordinates n, while the time will be denoted 
by uq. The link variables Up are matrices in color and have spatial matrix 
elements 

( u o(n )) mn = U (m,m )5 mtn 

( U j( m o)) mjn = ^(m,m )<Wj,n. (2) 

The quark field ip carries Dirac, color and flavor indices, denoted altogether 
by a, and space-time label n. We will use with the same convention Grass- 
mann variables x's and y's and quark creation and annihilation operators 
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x + } x,y + ,y. But unlike Grassmann variables, the creation and annihilation 
operators do not carry a temporal index. Thus according to our convention 

ip(n )ip(no) = 53Vn,o( 7l o)V'n 1 a(no) 

n.Q 

ry ' ry \ ryt ' nn 

n,a 

X + X(n ) = J2%ta X nA n o) (3) 

and the quark action is 

mo, no 

= J2 ^m,a( m o)Qma,nA m 0^ )^nA n o)- ( 4 ) 

m,a,n,/3 

Let us come back to the evaluation of the quark determinant in a given 
baryonic sector. 

The first step is to write the unconstrained determinant as a trace in Fock 
space 

det Q = Trf. (5) 

T is the transfer matrix whose definition will be given in the next Section. 
The second step is to impose the restriction to a sector with baryon number 
tib by inserting in the trace the appropriate projection operator P riB 

detQ\ nB =Tr(TP riB ). (6) 

To rewrite the above expression as a determinant we choose coherent states 
as a basis 

\x >= | exp(— xx + ) >, \y >= \ exp(—yy + ) > . (7) 
In such a basis the trace takes again the form of a Berezin integral 

detQ\ nB = J [dx + dxdy + dy] exp x + x — y + y^j < x,y\TP nB \ — x, —y > 

(8) 

which will be expressed in terms of the determinant of a modified quark 
matrix. 
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So far for the strategy. Now we show how we get a series expansion, what 
is its the physical interpretation and what is the final result. 
The kernel < x, y\TP nB \ — x, —y > has the integral form 

< x, y\TP„ B | — x, —y > = J [dz + dz dw + dw] exp z + z — w + w^j ■ 

< x, y\T\z, w >< z, w\P nB I — x, —y > . (9) 

The kernel < x,y\T\z,w > will be given in the next Section, while that 
of P ng is immediately calculated 

OO 1 

<z,w\P ng \-x,-y> = E(- 1 ) ns 77 T \ i n2 ' 

{{n B + r)\r\y 

< {yw + ) r {xz + y nB+r \xx + ) nB+r {yy + y > (10) 

Since x + , y + are creation operators of quarks and antiquarks respectively, we 
see that the r-th term of this series gives the gauge-invariant contribution of 
n B valence quarks plus r quark-antiquark pairs. 

Needless to say, forn B = 0, det Q\ nB does not reduce to the unconstrained 
determinant: Indeed also baryonic states are present in the unconstrained 
determinant, while they are not in detQ\ nB= o- In QCD at nonvanishing 
temperature it makes a difference whether we impose or not the condition 
riB = 0. In view of the relatively low value of the critical temperature with 
respect to the nucleon mass, however, we do not expect significant effects 
from this restriction unless we go to exceedingly high tempreatures. 

Let us now proceed to derive our final result. By evaluating the vacuum 
expectation values appearing in the last equation we express the kernel of 
the projection operator in terms of Grassmann variables only 

OO 1 

< z,w\P nB \-x,-y >= E(- 1 )" B , T \"j . (z + xy B+r (w + yy. (11) 

^ (n B + r)\r\ 

To evaluate the integral of Eq.(||) we rewrite the above equation in exponen- 
tial form 

°° 1 Q n B+r Qr 

< Z,w\P nB \ X ,y >= E + r) , r , dj n B +r gjr ^M~3^ + X - J 2 W + y)\ jl=J2=0 .{12) 

The integrals of Eqs.(@),(j|) are gaussian and we get the constrained de- 
terminant in terms of the determinant of a modified quark matrix 

°° 1 Q n B+r fir 

<•* oi« = £ (ns+ r)Mafi~m {Q + SQ) ^ (13) 
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The variation of the quark matrix is 

5Q = K[( Jl -l)(l + 7o) C /(iVo-l)t/o(iVo-l)r o (+) + 

(j 2 - 1)(1 - 7o) d(N ) Tt ] U+(N - 1) 

+ (hj 2 - 1) d(N ) C{N ) (1 + 7 o)] , (14) 

where 

Ui)^ _ r c 

I ) m „-|-i n „0 m n , 

V / m,n 

(d(N )) mn = 5 m0)No 5 mQ)no 5 m ^ (15) 

K is the hopping parameter and the matrix C is given in the next Section. 
The particular time N appearing depends on the fact that we inserted the 
projection operator at the latest time, but the value of detQ\ nB does not 
depend on this arbitrary choice, as it will be clear from the derivation. 



3 Evaluation of the quark matrix in a given 
baryonic sector 

In the first part of this Section we report the definition of the transfer matrix 
T and the proof of Eq. (|). These results have been derived by Luscher ||, 
and we will follow his paper even in the notation with minor changes. Then 
we will evaluate the integral (^) getting the quark matrix in a given baryonic 
sector. 

To start with we report the explicit expression of the quark matrix 

Q = K[(1 + 7o)^oT (+) + (1 - lv)T^U+ + 2C] - B (16) 

where 



Z 3=1 

B = t-Kj2(u j + U+) (17) 



and the lattice spacing has been set equal to 1. The modifications necessary 
at nonzero temperature are obvious. The transfer matrix is defined in terms 
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of the operator 

f F {n ) = exp{2Kx + B{n )-^c{n )B(n )-^y- 

■exp(-x + M(n )x-y + M(n )y^) (18) 
according to the ordered product 

T = J I] ( T ^K)) T F (n + l). (19) 

n =-N 

J is the jacobian of a transformation which will be defined later, and two 
new matrices have made their appearance [] 

M = -]n((2K)*UfB-*) 
c = l't(U s -Uf)i* j . (20) 

Z 3=1 

The (7j are the Pauli matrices and the operator Tp(no) depends on the time 
n only throu the dependence on it of the gauge fields. 

We then use the following equations. Firstly, for arbitrary matrices M 
and N 

< x\ expx + M x\x' >= exp (x + expMx'j , (21) 



< x\ exp(x + M x) exp(x + iV x)\x' >= exp (x + e M e N x') . (22) 

Secondly, if B = B(x + ) and C = C(x) are operators which depend on x + ,x 
only, for any operator A 

< x\B ■ A ■ C\x' >= B{x + )A{x + , x')C{x'). (23) 

We thus get the kernel 

< x(n ),y(n )\f^(n )f F (n + l)\x(n + l),y{n + 1) >= 
exp (2Kx(n )B(n y^c(n ) B(n )~^y(n )) 

■ exp (a; + (no)e M(no) e M(no+1) a;(n + 1) - y + (n )e M(no) e M{na+1) y(n + 1) 
• exp (2Ky{n + l)B(n + l)~h(n + 1) B(n + l) _ ^(n + 1)) . (24) 



1 Luscher assumes the gauge Uo = 1 to prove reflection positivity. We are not concerned 
here with this most important issue and do not fix the gauge. 
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Now to to go back from the Grassmann variables x 1 and y 1 to the quark 
field ip we only need the relations 

i> = (x + i/) 7o 5-i (25) 

This transformation cancels out the jacobian J . Inserting the identity be- 
tween the factors in Eq. (|I9"D , using the equality 

W)„ = X t X n + VnVn, (26) 

and collecting all the exponents we get Eq. (|5]). At this point two observa- 
tions are in order. The first is that the interpretation of annihilation 
operators of quarks and antiquarks follows from the action of the charge 
conjugation on the quark field 

U = -Ctf, (27) 
where, with Liischer's convention for the 7-matrices, 

C = 7o72- (28) 
In fact these tranformations imply 

x = —ioiV 

y' = ia 2 x. (29) 

The second observation is that all the above construction requires the Wilson 
parameter r to be equal to 1. 

There remains the evaluation of Eq. (|]). We notice that 

< x,y\f\z,w >= rexp (x + (iVo - l)Dz + y + (N - l)D*w + wEz) , (30) 

where r does not depend on z, w, the star means complex conjugation and 

D = B-*(N -1)U (N -1)B-*{N ) 

E = B-^(N )c(N )B- l i(N ). (31) 
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We must therefore perform the gaussian integral 

< x,y\TP nB \ — x, —y >= J [dz + dz dw + dw]r exp z + z — w + w^j ■ 

exp (-jiz + x-j 2 w + y + x + (N - l)Dz + y + (N - l)D*w + w E zj 
= rexp (-jiX + (iVo - l)D x - j 2 y + (N - l)D*y + Jlj2 y E x) . (32) 

Comparing the kernel of the unconstrained transfer matrix with the above 
and taking into account the antiperiodic boundary conditions of the quark 
field in euclidean time we get the expression of Eq. fll4] ) for 5 Q. 

4 Conclusion 

In principle by the present approach one can perform investigations of QCD 
in any definite baryonic sector at any temperature, starting from baryon 
number zero up to high baryonic densities. In practice, apart from the lowest 
baryonic numbers, we are confined to small physical volumes. With the 
volumes actually accessible in numerical simulations, perhaps this is not a 
too severe limitation, and we can hope to get reasonable results with the first 
few terms at least for not too high density and temperature. 

One could think of different applications of our expansion by further spec- 
ifying or changing the projection operator. In the case of baryonic number 
1, for instance, by evaluating the expectation value of the total angular mo- 
mentum one can disentangle the valence and cloud contributions to the 
spin of the nucleon. Remaining to the lowest baryonic numbers, one might 
try to attack the study of the deuteron. 

Another obvious field where our approach is potentially relevant is the 
QCD phase diagram. In particular one could investigate the various super- 
conductive phases with spontaneous breaking of the color symmetry (see for 
instance M) by replacing the projection operator by the trial state of the 
BCS theory. 
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